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Abstract 

We have searched for an exotic decay of orthopositronium into a single 
photon and a short-lived neutral boson in the hitherto unexplored mass 
region above 900 keV / c 2 , by noting that this decay is one of few remaining 
candidates which could explain the discrepancy of the orthopositronium 
decay-rate. A high-resolution measurement of the associated photon en- 
ergy spectrum was carried out with a germanium detector to search for 
a sharp peak from this two-body decay. Our negative result provides the 
upper-limits of 2.0 x 10 -4 on the branching ratio of such a decay in the mass 
region from 847 to 1013 keV/c 2 , and excludes the possibility of this decay 
mode explaining the discrepancy in the orthopositronium decay-rate. 
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Over the last seven years, three precision measurements of the orthopositron- 
ium (o— Ps) decay-rate were performed, which report the decay-rate faster than 
the QED prediction [|TJ by 1400, 1300 and 920 ppm, respectively, corresponding 
to 6.7, 6.3 and 4.3 standard deviations 0. In attempts to explain the origin of 
this discrepancy, various exotic o— Ps decay modes have been investigated with- 
out obtaining any evidence so far. The decay mode into invisible final states 
was ruled out at 3 ppm M. The decay into two and four photons, forbidden 
by space-rotational invariance and by the charge-conjugation invariance of QED, 
was excluded at 233 and 8 ppm, respectively The decay into two photons 

and a long-lived, weakly-interacting vector boson was excluded at 10 ppm ||. 
The decay into a photon and a long-lived, weakly-interacting particle was also 
ruled out at 10 ppm for a mass less than 1010 keV/ c 2 [^], while that into a photon 
and a short-lived neutral boson with a mass less than 900 keV/ c 2 was excluded 
at 400 ppm || [], |HJ| . One of very few remaining candidates for the source of 



the discrepancy is the decay into a single photon and a short-lived neutral boson 
X° in the hitherto unexplored mass range above 900 keV/c 2 . We report in this 
letter a dedicated search for such a decay. 

A schematic view of the experimental setup is shown in Fig. 1. A 22 Na 
positron source of intensity 0.7 /iCi (2 mm spot diameter) is sandwiched between 
a plastic scintillator (NE104) of thickness 100 fim and a mylar sheet of thickness 
25 /im, each with a diameter of 12 mm. The source and a target of silica aerogel 
(dimensions 10 x 10 x 5mm 3 , density 0.1 gem -3 ) are held at the center of a pipe 
(18 mm diameter) made of 20/xm thick aluminized mylar. The pipe is filled with 
N 2 gas at 1 atmosphere. Most of the positrons emitted toward the target pass 
through the scintillator, giving light pulses to two photomultipliers (Hamamatsu 
H-3165-PV). When stopped in the silica aerogel, they form positronium-atoms 
for about 20% of cases. The total efficiency of o— Ps formation in our set-up was 
measured to be 3.0% per /5 + decay. 

The energy of the photon from the decay of o— Ps is measured by a planar 
high-purity germanium detector (diameter 16 mm, thickness 10 mm, Ortec GLP- 
16195/10-P). The germanium detector has a thin Be window of thickness 0.13 
mm, which results to a high efficiency to low-energy photons down to 5 keV. The 
energy resolution and the absolute peak efficiencies of the germanium detector are 
determined as a function of the photon-energy by using the line 7 peaks from var- 
ious sources of known strength, i.e., 22 Na, 57 Co, 133 Ba, 152 Eu, 210 Pb and 226 Ra, 
placed at the source position. The energy resolution obtained is 270, 525 and 
828 eV FWHM at 14.4 , 136.5 and 356.0 keV, respectively. 

The trigger and the data-acquisition system is arranged as follows. The pulses 
from the two photomultipliers are discriminated individually and the coincidences 
between them provide the start signals to the time-to-digital converter (TDC ; 
LeCroy 2228A), as well as the trigger signals to the CAMAC system. One out- 
put from the preamplifier of the germanium detector is fed through a fast-filter 
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amplifier (Ortec 474) with an integration-time of 100 nsec, into the peak-hold 
analog-to-digital converter (ADC-1 ; Hohshin C008), and also into a discrimina- 
tor whose output is used as the stop signal for the TDC. The other preamplifier 
output is amplified by a shaping amplifier (Ortec 572) with an integration-time 
of 3 usee. One output of the shaping amplifier is fed into the peak-hold analog- 
to-digital converter (ADC-2 ; Hohshin C011), to record the information of the 
low-energy region below 160 keV. The other output of the shaping amplifier is 
attenuated and fed into the other input of the ADC-2, to record the energy infor- 
mation in a wider range up to 600 keV. If the TDC stop signal is not generated 
within 5 //sec after the trigger, the data-readout is terminated and the system is 
cleared to accept the next trigger. 

The data were collected from 10 runs for a total of 2.26 x 10 6 sec. Throughout 
the data-taking period, the room temperature was controlled to ±0.4°C. The 
energy calibrations and the measurements of the energy resolution were carried 
out at intervals of 7 days. In addition, the position and the width of the 511 keV 
annihilation line were monitored by offline analysis of the data, which showed a 
stable peak position within 0.03 keV and a constant energy-resolution throughout 
the data-taking period. 

The offline selection and calibration process is as follows. In order to eliminate 
the pile-up effects, we demand the consistency between two ADC values with 
different integration-time : If the ratio of the two values deviates from the central 
value by more than 6%, the event is rejected as being affected by the pile-up. 
This selection has an efficiency of 92%. The time walk of the discriminator 
output is corrected by using the pulse height information from the ADC-2. The 
nonlinearlity of the TDC is corrected by using a flat spectrum between random 
triggers and clock signals. 

Figure 2 shows the time spectrum between the scintillator and germanium 
signals. A sharp peak of the prompt annihilations is followed by the exponential 
decay of o— Ps and subsequently by the constant accidentals. The time spectrum 
fits well to an exponential function plus a constant term. A decay time of 134.2 
±0.4 nsec is consistently obtained for various time spans if the fitting starts 
more than 150 nsec after the prompt peak. To obtain a pure sample of the o— Ps 
decay, we select the events in the time-window between 150 and 400 nsec after 
the prompt annihilations. The accidental contribution is measured by using the 
events in the time-window between 2000 and 3200 nsec after the prompt peak. 

Figures 3a and 3b show the measured energy-spectra from the o— Ps decay 
(o— Ps spectrum), thus selected and the accidental contribution being subtracted. 
The two-body decay into a photon and X° (7X decay) would appear in this 
0— Ps spectrum as a narrow peak on a smooth background. The peak position 
k p is related to the X° mass, m x o, by 



(m x o/2m e f 



3 



where m e is the electron mass. 



Peaks in the o— Ps spectrum have been searched for by scanning k p with a 
step of 0.088 keV from 10 to 160 keV, corresponding to the mxo region from 847 
to 1013 keV/ c 2 : At each k p , the o— Ps spectrum was fitted with the function 

S(k) + n(kp)P(k) , 

where S(k) and P(k) represent the smooth background-spectrum and the nor- 
malized peak-function, respectively. The fitting parameter n(k p ) corresponds to 
the number of photons from the 7A decay events under the peak at k p . The 
background-spectrum S(k) is obtained at each k p by fitting the o— Ps spectrum 
(excluding the region within 2.5 times the FWHM of the peak position) with 
a polynomial function of up to seventh order. Since the intrinsic width of X° 
is expected to be extremely small, the observable width of the peak would be 
dominated by the energy-resolution of the detector and by the thermal motion 
of o— Ps. The thermal kinetic energy of o— Ps has been measured in an aerogel 
similar to the one used in this experiment and is known to decrease to about 0.04 
eV within 60 nsec after formation [[D], 12]. The Doppler-broadening due to this 



thermal kinetic energy would contribute negligibly to the peak-width. Therefore 
the shape of the peak is assumed to be Gaussian, of a width determined by the 
energy-resolution of the germanium detector. 

No statistically significant peak beyond 4 standard deviations was detected 
in the o— Ps spectrum, resulting to upper limits on n(k p ). The limits on n(k p ) 
can be converted into the corresponding limits on the branching ratio of the 7 AT 
decay, Br(k p ), by using the formula 

Br(k v )- 



^(^p) -^ops 



where e(k p ) is the absolute peak efficiency for a 7-ray of energy k 



N ops represents the total number of the 0— Ps decays occurring within the 
150 - 400 nsec time-window, and is determined by fitting the "compton-free" 
region of the 0— Ps spectrum (between 400 to 490 keV, where the three-photon 



decays dominate) with the theoretical spectrum JR] of the three-photon decay 
folded with the energy- dependent efficiency It should be noted that Br(k p ) 
thus obtained has a small systematic error, being dependent only on the relative 
values of the efficiencies, not on their absolute values. 

The resulting upper-limits at a 90% confidence level on the branching ratio of 
7A decay are shown in Fig. 4 , as a function of mxo ■ As indicated in this figure, 
new upper limits on the 7A branching ratio in the mass range of X° from 900 
to 1013 keV/ c 2 have been obtained. Since our upper limits are at least a factor 
of 7 below the level which would explain the reported discrepancy of the 0— Ps 
lifetime, we conclude that the 0— Ps decay mode of 7A in the mass region below 
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1013 keV/c 2 can not explain the reported discrepancy of the orthopositronium 
decay-rate. 

We wish to thank Dr. S. Asai for suggestions and discussions, Dr. X. Zhang 
for kindly lending us the germanium detector, Mr. Y. Nakamura, Mr. N. 
Kimura and Mr. M. Matsui, of Japan Radioisotope Association, for preparing 
the positron source. The analysis was performed on the RS/6000 workstations 
supplied for partnership program between ICEPP and IBM Japan Ltd. 
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Figure captions 



Fig. 1. A schematic view of the experimental setup. Circular inset: a 
magnified view of the source region. 

Fig. 2. The time spectrum between the scintillator and germanium detector 
signals. 

Fig. 3. The o— Ps spectra after subtracting the accidental contribution, (a) 
Below 160 keV. (b) Up to 600 keV. 

Fig. 4. The resultant upper-limits at 90% CL. on the branching ratio of 7X 
decay in comparison with the existing limits ||. 
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This figure "figl-l.png" is available in "png" format from: 
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Over the last seven years, three precision measurements of the orthopositron- 
ium (o— Ps) decay-rate were performed, which report the decay-rate faster than 
the QED prediction [1] by 1400, 1300 and 920 ppm, respectively, corresponding 
to 6.7, 6.3 and 4.3 standard deviations [2]. In attempts to explain the origin of 
this discrepancy, various exotic o— Ps decay modes have been investigated with- 
out obtaining any evidence so far. The decay mode into invisible final states 
was ruled out at 3 ppm [3]. The decay into two and four photons, forbidden 
by space-rotational invariance and by the charge- conjugation invariance of QED, 
was excluded at 233 and 8 ppm, respectively [4, 5]. The decay into two photons 
and a long-lived, weakly-interacting vector boson was excluded at 10 ppm [6]. 
The decay into a photon and a long-lived, weakly-interacting particle was also 
ruled out at 10 ppm for a mass less than 1010 keV/c 2 [7], while that into a photon 
and a short-lived neutral boson with a mass less than 900 keV/c 2 was excluded 
at 400 ppm [8, 9, 10]. One of very few remaining candidates for the source of 
the discrepancy is the decay into a single photon and a short-lived neutral boson 
X° in the hitherto unexplored mass range above 900 keV/c 2 . We report in this 
letter a dedicated search for such a decay. 

A schematic view of the experimental setup is shown in Fig. 1. A 22 Na 
positron source of intensity 0.7 fid (2 mm spot diameter) is sandwiched between 
a plastic scintillator (NE104) of thickness 100 fim and a mylar sheet of thickness 
25 fim, each with a diameter of 12 mm. The source and a target of silica aerogel 
(dimensions 10 X 10 X 5 mm 3 , density 0.1 g cm -3 ) are held at the center of a pipe 
(18 mm diameter) made of 20/xm thick aluminized mylar. The pipe is filled with 
N 2 gas at 1 atmosphere. Most of the positrons emitted toward the target pass 
through the scintillator, giving light pulses to two photomultipliers (Hamamatsu 
H-3165-PV). When stopped in the silica aerogel, they form positronium-atoms 
for about 20% of cases. The total efficiency of o— Ps formation in our set-up was 
measured to be 3.0% per /3 + decay. 

The energy of the photon from the decay of o— Ps is measured by a planar 
high-purity germanium detector (diameter 16 mm, thickness 10 mm, Ortec GLP- 
16195/10-P). The germanium detector has a thin Be window of thickness 0.13 
mm, which results to a high efficiency to low-energy photons down to 5 keV. The 
energy resolution and the absolute peak efficiencies of the germanium detector are 
determined as a function of the photon-energy by using the line 7 peaks from var- 
ious sources of known strength, i.e., 22 Na, 57 Co, 133 Ba, 152 Eu, 210 Pb and 226 Ra, 
placed at the source position. The energy resolution obtained is 270, 525 and 
828 eV FWHM at 14.4 , 136.5 and 356.0 keV, respectively. 

The trigger and the data-acquisition system is arranged as follows. The pulses 
from the two photomultipliers are discriminated individually and the coincidences 
between them provide the start signals to the time-to-digital converter (TDC ; 
LeCroy 2228A), as well as the trigger signals to the CAMAC system. One out- 
put from the preamplifier of the germanium detector is fed through a fast-filter 
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amplifier (Ortec 474) with an integration-time of 100 nsec, into the peak-hold 
analog-to-digital converter (ADC-1 ; Hohshin C008), and also into a discrimina- 
tor whose output is used as the stop signal for the TDC. The other preamplifier 
output is amplified by a shaping amplifier (Ortec 572) with an integration-time 
of 3 fisec. One output of the shaping amplifier is fed into the peak-hold analog- 
to-digital converter (ADC-2 ; Hohshin C011), to record the information of the 
low-energy region below 160 keV. The other output of the shaping amplifier is 
attenuated and fed into the other input of the ADC-2, to record the energy infor- 
mation in a wider range up to 600 keV. If the TDC stop signal is not generated 
within 5 fisec after the trigger, the data-readout is terminated and the system is 
cleared to accept the next trigger. 

The data were collected from 10 runs for a total of 2.26 X 10 6 sec. Throughout 
the data-taking period, the room temperature was controlled to ±0.4°C. The 
energy calibrations and the measurements of the energy resolution were carried 
out at intervals of 7 days. In addition, the position and the width of the 511 keV 
annihilation line were monitored by offline analysis of the data, which showed a 
stable peak position within 0.03 keV and a constant energy-resolution throughout 
the data-taking period. 

The offline selection and calibration process is as follows. In order to eliminate 
the pile-up effects, we demand the consistency between two ADC values with 
different integration-time : If the ratio of the two values deviates from the central 
value by more than 6%, the event is rejected as being affected by the pile-up. 
This selection has an efficiency of 92%. The time walk of the discriminator 
output is corrected by using the pulse height information from the ADC-2. The 
nonlinearlity of the TDC is corrected by using a flat spectrum between random 
triggers and clock signals. 

Figure 2 shows the time spectrum between the scintillator and germanium 
signals. A sharp peak of the prompt annihilations is followed by the exponential 
decay of o— Ps and subsequently by the constant accidentals. The time spectrum 
fits well to an exponential function plus a constant term. A decay time of 134.2 
± 0.4 nsec is consistently obtained for various time spans if the fitting starts 
more than 150 nsec after the prompt peak. To obtain a pure sample of the o— Ps 
decay, we select the events in the time-window between 150 and 400 nsec after 
the prompt annihilations. The accidental contribution is measured by using the 
events in the time-window between 2000 and 3200 nsec after the prompt peak. 

Figures 3a and 3b show the measured energy-spectra from the o— Ps decay 
(o— Ps spectrum), thus selected and the accidental contribution being subtracted. 
The two-body decay into a photon and X° ("fX° decay) would appear in this 
o— Ps spectrum as a narrow peak on a smooth background. The peak position 
k p is related to the X° mass, m_yo, by 

k p = m e 1 - (m x o /2m e ) 2 \ , 
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where m e is the electron mass. 



Peaks in the o— Ps spectrum have been searched for by scanning k p with a 
step of 0.088 keV from 10 to 160 keV, corresponding to the rn x ° region from 847 
to 1013 keV/c 2 : At each k p , the o— Ps spectrum was fitted with the function 

S{k) + n(k p )P(k) , 

where S(k) and P(k) represent the smooth background-spectrum and the nor- 
malized peak- function, respectively. The fitting parameter n(k p ) corresponds to 
the number of photons from the "fX° decay events under the peak at k p . The 
background-spectrum S(k) is obtained at each k p by fitting the o— Ps spectrum 
(excluding the region within 2.5 times the FWHM of the peak position) with 
a polynomial function of up to seventh order. Since the intrinsic width of X° 
is expected to be extremely small, the observable width of the peak would be 
dominated by the energy-resolution of the detector and by the thermal motion 
of o— Ps. The thermal kinetic energy of o— Ps has been measured in an aerogel 
similar to the one used in this experiment and is known to decrease to about 0.04 
eV within 60 nsec after formation [11, 12]. The Doppler-broadening due to this 
thermal kinetic energy would contribute negligibly to the peak-width. Therefore 
the shape of the peak is assumed to be Gaussian, of a width determined by the 
energy-resolution of the germanium detector. 

No statistically significant peak beyond 4 standard deviations was detected 
in the o— Ps spectrum, resulting to upper limits on n(k p ). The limits on n(k p ) 
can be converted into the corresponding limits on the branching ratio of the "fX° 
decay, Br(k p ), by using the formula 



t{kp) A/" ops 

where e(k p ) is the absolute peak efficiency for a 7-ray of energy k p . 

N ops represents the total number of the o— Ps decays occurring within the 
150 - 400 nsec time-window, and is determined by fitting the "compton-free" 
region of the o— Ps spectrum (between 400 to 490 keV, where the three-photon 
decays dominate) with the theoretical spectrum [13] of the three-photon decay 
folded with the energy-dependent efficiency. It should be noted that Br(k p ) 
thus obtained has a small systematic error, being dependent only on the relative 
values of the efficiencies, not on their absolute values. 

The resulting upper-limits at a 90% confidence level on the branching ratio of 
"fX° decay are shown in Fig. 4 , as a function of mx°- As indicated in this figure, 
new upper limits on the "fX° branching ratio in the mass range of X° from 900 
to 1013 keV/c 2 have been obtained. Since our upper limits are at least a factor 
of 7 below the level which would explain the reported discrepancy of the o— Ps 
lifetime, we conclude that the o— Ps decay mode of "fX° in the mass region below 
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1013 keV/c 2 can not explain the reported discrepancy of the orthopositronium 
decay-rate. 

We wish to thank Dr. S. Asai for suggestions and discussions, Dr. X. Zhang 
for kindly lending us the germanium detector, Mr. Y. Nakamura, Mr. N. 
Kimura and Mr. M. Matsui, of Japan Radioisotope Association, for preparing 
the positron source. The analysis was performed on the RS/6000 workstations 
supplied for partnership program between ICEPP and IBM Japan Ltd. 
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Figure captions 



Fig. 1. A schematic view of the experimental setup. Circular inset: a 
magnified view of the source region. 

Fig. 2. The time spectrum between the scintillator and germanium detector 
signals. 

Fig. 3. The o— Ps spectra after subtracting the accidental contribution, (a) 
Below 160 keV. (b) Up to 600 keV. 

Fig. 4. The resultant upper-limits at 90% C.L. on the branching ratio of "fX° 
decay in comparison with the existing limits [8]. 
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